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Abstract 

As a contribution to the study of Hartree-Fock theory we prove rigor- 
ously that the Hartree-Fock approximation to the ground state of the d- 
dimensional Hubbard model leads to saturated ferromagnetism when the 
particle density (more precisely, the chemical potential fj,) is small and the 
coupling constant U is large, but finite. This ferromagnetism contradicts the 
known fact that there is no magnetization at low density, for any U, and thus 
shows that HF theory is wrong in this case. As in the usual Hartree-Fock 
theory we restrict attention to Slater determinants that are eigenvectors of 
the z-component of the total spin, S 2 = Yl,x n x,\ ~ n x,l> an< ^ we ^ nc ^ tnat 
the choice 2§ z = N = particle number gives the lowest energy at fixed 
< \i < Ad. 
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1 Introduction 

The (one-band) Hubbard model has become a standard model for correlated elec- 
trons in condensed matter physics since it is, perhaps, the simplest possible model 
of itinerant interacting electrons. In spite of its simplicity, its zero temperature 
phase diagram is rich with different magnetic phases such as paramagnetic, fer- 
romagnetic, and antiferromagnetic phases, depending on the details of the hop- 
ping amplitudes, the (relative) coupling constant U/t and the filling parameter 
i/ = JV/(2|A|). 

As the Hubbard model is a many-body fermion model, the computation of 
its ground state for large lattices is a difficult, if not impossible, task, except in 
one-dimension HI El. Thus various schemes have been developed during the past 
decades to derive an approximate ground state and then to study its magnetic phase 
diagram. 

In the present paper, we consider the Hartree-Fock approximation of the (re- 
pulsive, one-band, nearest-neighbor-hopping) Hubbard model with the intention 
of studying the validity of the Hartree-Fock approximation. We require the Slater 
determinants entering the Hartree-Fock energy functional to be eigenfunctions of 
the operator E> z : = ^LeA^^t ~ n x,i\ °f tota l s P m m tne -z-direction, and for this 
reason we refer to the model as the HFz approximation. Our requirement means 
that each orbital has the form ip(x) (g> |f) or (p(x) <8> ||). This is a restriction in 
the sense that general orbitals are of the form ip(x, a), in which the spin direc- 
tion depends on position. No other restriction is imposed on the variational states; 
in particular, no assumption about translation invariance is made a priori. For 
the HFz model, at small chemical potential and for sufficiently strong repulsion, 
we give a mathematical proof of saturated ferromagnetism in the Hartree-Fock 
ground state. That is, the HF ground state has maximal total spin and maximal 
ferromagnetic long-range spatial order. The smallness of the chemical potential 
and the large strength of the repulsion also insure that the HF ground state density 
is strictly below half- filling. 

Before we come to a detailed description of our result and its proof, we discuss 
it in comparison to other works. 

The appearance of ferromagnetic behaviour has been anticipated in many stud- 
ies of the Hubbard model and approximations thereof. Among these are (re- 
stricted) Hartree-Fock approximations Q, DMFT models in the limit of infinite 
spatial dimension |H1|5J|5J|7|, exact diagonalizations on small lattices (5), varia- 
tional calculations [9| and studies at low filling ifTUl . These studies support the 
conjecture that, for large coupling U/t ^> 1 and away from half-filling, v ^ 1/2, 



BUT, Version of 24-M-2006 



3 



the ground state of the Hubbard model is ferromagnetic. Ferromagnetism has 
been established for the (full) Hubbard model in case the dispersion relation leads 
to a very high density of states around the Fermi energy ifTTl fT2l [T3ll and in case 
of next-nearest-neigbor hopping [fT4l FTSl . 

As said before, the main purpose of the present paper is to prove ferromagnetic 
behaviour with mathematical rigor. None of the papers |3J |5J [TJ [HI cited above 
match the standards of a mathematical proof: The orbitals in the Hartree-Fock 
approximation are a priori assumed to be composed of only few Fourier modes; 
the error terms when taking the limit of infinite spatial dimension in DMFT are 
not under control; exact diagonalizations are restricted to very small lattices and 
the implication of these to the thermodynamic limit remains unclear. The work 
by Mielke and Tasaki lfTTl[T^l[T3l is mathematically rigrous, but the assumptions 
made therein about the lattice structure are rather special. On the other hand, 
by adding next-nearest-neigbor hopping (two-band Hubbard model), Tasaki lfl4l 
15 1 has found a Hubbard model that displays ferromagnetism in all dimensions. 
Tasaki also reviews rigorous results on ferromagnetism in the Hubbard model in 

m. 

While the prediction of ferromagnetism in the Hubbard model and approxi- 
mations thereof is supported by the above studies, we also know that HF theory 
predicts anrf-ferromagnetism (in the sense that the total spin is zero) at higher 
densities, notably at half-filling [fTTl. Furthermore, our proof shows saturated fer- 
romagnetism at low density and sufficiently large coupling in HF theory, even in 
one-dimension, but the actual ground state always has spin zero in one-dimension 
as long as there is only nearest-neighbor hopping (see CHI). 

Even more seriously, our conclusion is opposite to what actually occurs in the 
Hubbard model. Namely, at very low density (and independent of the value of 
U > 0), there is no magnetization in the ground state of this model. In the ground 
state S 2 is close to zero and converges to zero, as the particle density tends to zero. 
This has been pointed out in |[T9lfT6l . based on arguments similar to the following 
transcription to lattice systems of the recent work ll20ll . 

In this paper [ 20 1 it was shown that fermions in the 3-dimensional continuum 
M 3 (instead of the lattice Z 3 ), and with a repulsive two-body potential, have a 
ground state energy density, e, given by 

e(phPl) = ^5( 67r2 ) 2/S ( p T /3 + P| /3 ) + ^ 87rc Wl+ m g her order in (p T ,pJ , 

(1.1) 

where p-f, pj. are the densities of the 'spin-up' and the 'spin-down' fermions and a 
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is the scattering length of the two-body potential. Because p 5 ^ 3 dominates p 2 for 
small p, it is clear from (11.11) that the minimum energy occurs approximately, if 
not exactly, when pj = pj = p/2. This answers the questions in ET1 problem 3]. 

To show that there is vanishing net magnetization as p — > one only needs an 
upper bound for e of the form dl.ll) . For the Hubbard model (where the two-body 
potential is a positive delta-function, or even a hard core) this can conveniently be 
done by a variational wave-function of the form \P = F^ , where \l/ is a Slater 
determinant, and F is the projection onto the states with no double occupancy 
- in imitation of lfT9lfT6ll20ll . We omit the details, but we draw attention to the 
fact that F^ is not a Slater determinant, reflecting the more complex structure 
of correlations in the actual ground state of the Hubbard model. The proof of an 
analog of dl.lt with precise constants is a more complicated matter which is now 
under investigation, but it is not needed for the present discussion. 

Our setting is the usual (repulsive) Hubbard model with nearest-neighbor hop- 
ping on a d-dimensional cubic lattice A, with periodic boundary conditions and 
linear size L, which we assume to be an even integer. It is defined by the second 
quantized Hamiltonian 

H^U = Yl (~ A ^,y ~ l 15 ^) C l,a C y,a + ™*,T n *,l ■ ( L2 ) 

x,y£A, <t=T,J. x-GA 

We work at fixed chemical potential p instead of fixed particle number. The only 
slightly unusual notation is A x >y = T x>y — 2dS x>y for the matrix elements of the 
discrete Laplacian A on A, with T x>y := l[\x—y\i = 1] being the nearest-neighbor 
hopping matrix and 6 xy — l[x — y] the Kronecker-Delta. 

The operators c* a , c x a , and n XjCT := c* a c x a are the usual fermion creation, 
annihilation, and number operators, respectively, at site x E A and of spin a 6 
{t, I}, obeying the canonical anticommutation relations {c x a) c y T } = {c* , c* T } 
= 0, {c X (7 ,c* T } = 5 XiV 5 ajT , and c X (J |0) = 0, for all x,y,a,r. Here |0) is the 
vacuum vector in the usual Fock space JF A := JF^(C A <g> C 2 ) of spin-| fermions. 
The Hamiltonian H^u depends parametrically on the chemical potential p > 
and the coupling constant U > 0. 

Note that the usual hopping parameter t equals 1 here and that the discrete 
Laplacian A differs from the usual hopping matrix by the inclusion of the diagonal 
term, i.e., 2d times the identity matrix. This difference amounts to a convenient 
redefinition of the chemical potential p, so that p = corresponds precisely to 
zero filling since the hopping matrix —A > is a positive semi-definite matrix. 
Moreover, the boundedness < p < 4c? of p together with the assumption that 
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U ^ Ad insures that the corresponding electron density in the HF ground state is 
always at low filling, i.e., strictly below half-filling, < p < 1. 

Our definition of p is convenient because in this paper, we are concerned with 
the Hubbard model at low filling, and Our assumption of a bounded chemical 
potential < fi < 2d 

Apart from this, everything is standard. 

The Hamiltonian H^ u is a linear operator on the Fock space and the ground 
state energy E^j is its smallest eigenvalue, 

:= min{(^| HV) | * G T k , ||*|| = l} . (1.3) 

As the dimension dim(^A) = 2 dim(cAcg,c2) = 4 (L<i) < oo is finite, the determina- 
tion of E^ amounts to diagonalizing the finite-dimensional, selfadjoint matrix 
H^ ij. The fast growth of this dimension with the number L d of points in the lattice 
A, however, allows for an explicit diagonalization of H^ u by a modern computer 
only up to L = 4, in three spatial dimensions, d = 3. 

The Hartree-Fock (HF) approximation is an important method to reduce the 
high-dimensional many-particle problem given by the diagonalization of H^u to 
a low-dimensional, but nonlinear variational problem. It is defined by restricting 
the minimization in (11.31) to Slater determinants (pi A • • • A (fx, where {ipi}f =l C 
C A ® C 2 is an orthonormal family of N one-electron wave functions. The HF 
approximation to the Hubbard model was analyzed in ifTTl in the special situation 
when the number of electrons equals the number of lattice sites, N = | A|, which 
is usually referred to as half-filling. 

Note that a priori no other condition but orthonormality is imposed on the or- 
bitals {Lpi]f =l in the Slater determinants varied over in Hartree-Fock theory. This 
is sometimes stressed by calling it the unrestricted Hartree-Fock theory. Let us 
temporarily consider a general many-body Hamiltonian H which commutes with 
a certain symmetry operator S, i.e., [H, S] = 0. It is important to note that in 
this case, the HF ground state &hf, i-e-, the Slater determinant which minimizes 
the energy (Q>hf\H &hf), is not necessarily an eigenstate of S. Phrased differ- 
ently, unrestricted Hartree-Fock theory may (depending on the model) break the 
symmetry S. The following are examples that occur in physically relevant situa- 
tions: unrestricted HF ground states of atoms are, in general, not eigenfunctions 
of the angular momentum operator (because in unrestricted HF theory, all shells 
are filled [ 22 1) - even though the atomic Hamiltonian is rotationally invariant; the 
ground state in the BCS theory of superconductors (which is a variant of HF the- 
ory) is not an eigenfunction of the number operator - even though the BCS Hamil- 
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tonian preserves the particle number; a HF ground state for the Hubbard model 
with non-zero spin breaks the invariance of the Hubbard Hamiltonian under global 
spin rotations; charge density waves (CDW) and spin density waves (SDW) of the 
Hubbard model are translation invariant only by translation of an even number 
of lattice sites, breaking the (full) translation symmetry the Hubbard Hamiltonian 
Hn,u posesses. As it is impossible to predict a priori whether a symmetry of the 
Hamiltonian is preserved or not, we call all variations of (Q\H $) over Slater de- 
terminants $ which fulfill an additional constraint restricted Hartree-Fock theory. 

In this paper, we consider a restricted Hartree-Fock theory, which we term 
the HFz approximation. The further restriction imposed is that we minimize 
in (11.31 ) only over Slater determinants $ that are eigenfunctions of the operator 
S 2 := XLeA{ n ;c,T — n x ^} of total spin in the z-direction. One could rephrase our 
condition by saying that we do not allow for spiral spin density waves (SSDW; 
see, e.g., 0) in d 1 .31) . Once again, it is customary to employ this restriction in 
HF calculations without explicitly drawing attention to the fact that this is a re- 
striction. (In ifTTll mentioned above, however, we dealt with truly unrestricted HF 
theory.) 

More concretely, our HF wave functions have the form 

iV T N t 

$ = II c T^)n c I(^)iO}, (1.4) 

1=1 3=1 

where c| , (/) = J2 x &a f( x ) c z,t,J.' tne mte S ers are me particle numbers, and 
where the fi and gi are two families of orthonormal wave functions on the lattice 
A, i.e., (fi\fj) = (gi\gj) = 6 itj , with (f\g) : = ^ xeA f{x)g{x) denoting the usual 
hermitian scalar product for such functions. 

It is convenient to rephrase the HFz approximation in terms of one-particle 
density matrices, i.e., complex, self-adjoint A x A matrices whose eigenvalues lie 
between and 1 . To this end, we denote 

:= -A -/i (1.5) 

and observe that 

i=l j=l 

JV-j N l 

+ ^E(Ei/^)i 2 )(Ei^)i 2 )- < L6 ) 

xeA i=l j=l 
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Introducing the one-particle density matrices 7^ corresponding to $ by 

7t := ^2\fi)(M and 7| := ^2\9i)(9i\, (1-7) 
i=i i=i 

we observe that 74,4 = 7* , = 7* , are orthogonal projections of dimension jVj^ 
and that the energy expectation value of the Slater determinant $ is given by 

($\H$)=£^\ lhll ), where 

4!£ Z) (7t>7i) := Tr{K M ( 7T +7|)} + 17 X>(*) PiO«0 > d-8) 

and the diagonal matrix elements Pi,i(x) := ( r Yi,i) x ,x of 74^4 are the one-particle 
densities of the electron with spin up ("T") and spin down ("j"), respectively. 

The symbol "Tr" denotes the usual trace Tr{v4} = XLeA A x ,x of a complex 
A x A matrix A = (A X:y ) XjyeA with A x y £ C. That is, "Tr" is the trace over the 
states in C A of a single spinless particle on the lattice A. It does not include spin 
states, and it is not the trace over states in Fock space. 

Let us note that the particle numbers jVj^ are not determined ab initio. We are 
in the grand canonical ensemble, so they are determined by the condition that the 
total energy (11.81) is minimized. 

These observations motivate us to define the HFz energy by the following 
variational principle over projections: 

E%> := min {^74, 7|) | lU = 7*,4 = 7?,|} ■ ^ 

The two sets of orthogonal projections on C A over which we minimize in (11.91 ) is 
not really well-suited for a variational analysis. In particular, they are not convex. 
An observation in ll23ll . however, states that, because U > 0, we will obtain the 
same value for the minimum if we vary over the larger set of all one-particle 
density matrices, < 74,4 < 1, not only over projections. (Recall that a density 
matrix is a hermitean A x A matrix 7 whose eigenvalues lie between and 1, i.e., 
< 7 < 1, as a matrix inequality.) Our extended E^ff is then 

E$> = min{4^( 7T ,7|) | < 74,4 < 1} ■ d-10) 

The evaluation of E^ 1 ^ and the determination of those pairs (71,74.) of one- 
particle density matrices that minimize £^ is the objective of this paper. Our 

main result is that, for any < /x < Ad, the minimal value of Sjffi is attained for 
the saturated ferromagnet, provided U < 00 is sufficiently large. 
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Theorem 1.1 (Ferromagnetism). For any < /i < Ad, there is a finite length 
L#(/i) and a finite coupling constant U#(u) > 0, such that, for all even L > 
L#(fi) and all U > U#(p), the minimal HFz energy is given by the sum of the 
negative eigenvalues of —A — fi, 

rffi? = Tr{[-A- M ]_}. (1.11) 

If fi is not an eigenvalue of —A and if (7-f , 7^ ) is a minimizer of the HFzfunctional, 
i.e., < 7 Ta < 1, on<i^y(7 T ,7 i ) = E { ™j\ then 

either 7 T = 1[-A < p] , j l = (1.12) 
or 7T = 0, 7| = l[-A</i], (1.13) 

where 1[— A < //] is ?/ze spectral projection of —A onto (—00, /x). 

[With reference to Eq. (11.1 II) and elsewhere, note that in our notation, = 
min{X, 0} is negative, whereas elsewhere one often defines [Jf]_ to be positive, 
i.e., [X)_ := max{— X, 0}. If X is a self adjoint operator then [Jf]_ denotes the 
negative part of X and Tr[X]_ is the sum of the negative eigenvalues of X.] 

Theorem ll.ll is not really as complicated as it looks. It is stated in terms of a 
length L # and coupling constant U# in order to make it clear that the state of sat- 
urated ferromagnetism is obtained not only asymptotically in the thermodynamic 
limit and asymptotically as U — > 00, but it holds for all systems with large, finite 
interaction and sufficiently large size. 

Theorem 11.11 states that, for any value of the chemical potential \i E (0, Ad), 
the HFz variational principle yields a ferromagnetic minimizer, provided U and L 
are chosen sufficiently large (but still finite). A similar statement was proved in 
ifTTl Theorem 4.7] for U = 00 (which amounts to requiring (&\n x ^n^^f) = 0, 
on every lattice site x E A). 

At first sight, Theorem 11.11 seems to contradict another fact proved in [ 17 1 that 
the HF minimizer is antiferromagnetic at half-filling. But as the definition of the 
chemical potential ji in present paper differs from its definition in IfTTl by 2d + U, 
the parameter range of the present paper and of IfTTl never overlap and, hence, 
there is no contradiction. 

As just mentioned, the minimal HF energy and the minimal HFz energy agree 
in the half-filling case, as shown in IfTTl . We conjecture that this is also the case 
for the range of the chemical potential fx E (0, Ad) and sufficiently large U, but 
we do not know how to prove this conjecture. This is a topic for future research. 
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From Theorem ll.il we conclude that at small filling there is a phase transition 
(within the context of HFz theory) from paramagnetism for small U to saturated 
ferromagnetism for large U . This follows from continuity and the fact that when 
U = we can find the ground state explicitly and, as is well known, it has 5 = 
and is obtained from filling up the Fermi sea for both f and j states. 

If < fi < | then we can estimate L#((jl) and U#((i) in Theorem l 1 . 1 I more ex- 
plicitly. For the precise formulation of these estimates, we introduce the following 
constants, 

L*(/j) := 2 MM : = 24 {Ad) V" 2 , (1.14) 



,d 



k(jm) := l + 2 1n(2)(r 1 + l)+ln(4^- 1 ) 



42d+l e d (jd 

^(2+d)/2 



-2d 



(1.15) 



aM '= ¥W* (2nY {Adf (L16) 

r / \ . ( a 2 a k(u) ^ 

o* u, a) := mm < — , — — , > (1.17) 

! l(12d) 2 ' 3/i[4M,(/x) + l] d ' 2 J V 

TK \ r 2 ^ 24rf 2 ^ 

£/*(//, a) := max — -, — — -f, (1.18) 

where (S^ 1 ) = 2n d / 2 /T(d/2) is the measure of the unit sphere in Mr. 

Theorem 1.2. For any < /i < |, Theorem \l.l\ holds true with L#(p) := 
and U#(/jl) := as defined in A1.14\) . A1.16\) . and rti.iffl ). 

The explicit form of a*(/i), and £/*(/x, «*(//)), for a given < /i < |, 

in Theorem II .21 allows us to estimate the actual minimal size of L and U that 



guarantees saturated ferromagnetism. The distinction between /i < 1/2 and /i > 
1/2 is not a fundamental one. It is an artifact of the use in Lemma l3~6l of refs. ll2"4l 
and |[25|, whose methods favored this technical distinction. 



Acknowledgements: The authors are grateful to Alessandro Giuliani for very 
helpful discussions and comments about an earlier version of this paper. They 
also thank Manfred Salmhofer, Jiirg Frohlich, and Daniel Ueltschi for useful dis- 
cussions. MT thanks the german student exchange service DAAD for a generous 
stipend, which supported two thirds of his graduate studies. VB and MT grate- 
fully acknowledge financial support from grant no. HPRN-CT-2002-00277 of the 
European Union and grant no. Ba 1477/3-3 of the Deutsche Forschungsgemein- 
schaft. EL gratefully acknowledges support from the Alexander von Humboldt 



BLT, Version of 24-M-2006 



10 



Foundation of a fellowship, the U.S. National Science Foundation, grant no. PHY- 
0133984, and the hospitality of the Mathematics Departments of the University of 
Mainz and the Technical University of Berlin. The authors appreciate the careful 
and helpful work of a referee. 



BUT, Version of 24-M-2006 



11 



2 Proofs of Theorems QUI and 11.21 

This section contains the proofs of our main results, Theorems 11.11 and 1 1 . 21 with 
the aid of several lemmas which will be proved later in Section |3J Here is a brief 
outline of the strategy of the proof. 

• We first reduce the minimization of E^ff (7^, 7jj in (11.101) over two one-particle 
density matrices 7^ and j± to the minimization of an effective energy functional 

(7) which depends only one one-particle density matrix 7. It is given as a 

sum of two terms, £^u\ r y) = Tr-fliyy} + Tr{[i^ /t + U p]-}, where we recall that 
K^ = -A-p. 

• Given a trial one-particle density matrix 7 and a small number 5 > 2pU~ 1 , we 
introduce the corresponding particle density p{x) := ^ x>x and define the regions 
Vt := {x\p(x) < 5} and f2 c := {x\p(x) > 5} of low and high density onto which 
we project by Pq = J2 x en \ x ) ( x \ anc ^ = t — Pn, respectively. 

• We then use the fact that 7 is mostly localized in the high density region Vt c . 
This leads us to estimate the kinetic energy Tr{— APq'jPq} in f2 by zero and 
Tr{— APq 7PQ } in Q c by the kinetic energy of the free Fermi gas in fi c . The 
localization error is of order of a small constant times the volume \dfl\ of the 
boundary of fi. In Lemma | 3.1l we give the exact formulation of the bound which 
we use to estimate the term Tr{iiyy} in E^ (j). 

• For the analysis of the term Tr { [K^ + Up] _ } in Ejffi (7) , we use the fact that fl c 
is a classically forbidden region, because —p + Up > —p + U5 > p in tt c . So, as 
shown in Lemma l3~2l we can replace Tr{ [K^+Up] _ } by Tr{ [P n (K IM +Up)P n ] -}, 
up to localization errors of order of a small constant times \dQ\. 

• We then pick a (large, but fixed) number M > 1 and further split up the low 
density region £1 into the subset £li of those points in Q that are at most at distance 
2M away from the boundary dVl and the bulk Q 2 C Cl of points of distance 2M or 
more to <9S1 The contribution of f2i turns out to be negligible because fi x contains 
at most (AM + l) d |<9f2| points, and the density is low in Qi C £7. 

• The estimate of the region Q 2 3 x then uses the lower bound on the spatial 
density ^[K^ + Up < 0](x, x) of the projection onto the negative eigenvalues of 

+ Up (actually, p instead of p), which we derive in Lemma 13 .31 

• Adding up the estimates derived so far, we finally observe that £^u (j) is 
bounded below by Tr{[P n ^P n ]_} + Tr{[P£K„P£]_} - v \dn\ =:Y — r}\dQ\, 
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where r] > becomes small when U ^> 1 and 5 > is properly chosen. 
In Lemma | 3. 61 we reproduce the result from ll24l l25ll that Y can be estimated 
from below by Trjfi^].} + a|<9f2|, where a > depends only on p. In other 
words, the introduction of a domain wall at OVl drives up the energy by ct|<9f2|, 
which dominates 7/|<9fi|, provided rj is small. This establishes that 8^ - (7) > 
Tr-Qi^].} + (a — r))\dCl\, which implies the claim. 

To carry out the proof in detail, we start with the observation that the mini- 
mization over two one-particle density matrices in (11.101) can actually be reduced 
to the minimization over only one one-particle density matrix. To see this, we 
observe that 

$^Pt(X)P|(X> = Tr{p T7i }, (2.1) 
where p^ acts as a multiplication operator, (p^f) (x) := p^(x)f(x). Thus we have 



Tr{K^} + min (Tt {(K , + U p^) Tl }) 

0<7|<1 



(2.2) 



= o min i (Tr{/^ 7T } + Tr{[K, + Up,].}) . (2.3) 

(Recall that = —A — p.) In other words, we have done the minimization 
over 7 j in (12.21) by taking 7^ to be the projection onto the negative eigenspaces of 
+ Upi . Thus, as our minimization principle over only one 7, we obtain the 
following. 

^Effl = min{^ ) (7)|0< 7 <l}, (2-4) 
8^(j) := Tr{i^ 7 } + Tr{[K, + Up)„}, (2.5) 

where p(x) := 7^. From now on 7, with < 7 < 1, is an arbitrary, but 

fixed one-particle density matrix, for which we bound 8jj h Jj (7) from below. (An 
upper bound that agrees with Theorem ll.ll is readily obtained simply by choosing 
the variational function consisting of the unperturbed Fermi sea with all particles 
spin-up or all spin-down.) 

For the next step of the proof we introduce a small number 5 > 2pU~ 1 , whose 
precise value will be chosen in the final step of the proof. Given a one-particle 
density matrix < 7 < 1 with corresponding density p(x) := 7 X>X , we write the 
lattice A = U £l c as a union of two disjoing subsets of A in the following way. 

Q, ■= [ x e A I p(x) < 5} , (2.6) 
Vl c := {x e A I p(x) > 5} . (2.7) 



BLT, Version of 24-M-2006 



13 



These are the regions of low and high density, respectively. We define the bound- 
ary dVl of f2 by 

dtt := {x G n | disti(ar, fi c ) = 1} , (2.8) 

where disti(x, A) is the length of (number of bonds in) a shortest path joining 
x and some point in y G A. Another useful notion of distance which we shall 
use is distoo(x, A), which is defined by the condition that 2 distoo(a;, A) + 1 is the 
sidelength of the smallest cube centered at x that intersects A. When A is a single 
point y these distances are denoted by \x — y\i and \x — y]^. 

We define Pq, Pqc = P£, and Pqq to be the orthogonal projections onto f2, 
f2 c , and dil, respectively, where the projection onto an arbitrary set A C A is 
given by 

tut <*» 

We further set 

n(r) ■ I forxeW, 

P{X) - \ minj^j, p(x)}, for x Gil, (2JU) 

and observe that p(x) < p(x), for all x G A, which implies that 

^ z) ( 7 ) > Tr{K^} + Tr{[K„ + Up\_}. (2.11) 

For brevity, we define M := M*(//) := 12 (y) 2 and note that, by assumption, 
L obeys L > 2M. We further decompose f2 into two disjoint subsets and fi 2 
defined by 

fti := {x G fl | distoo(a;, Q c ) < 2M} , (2.12) 
n 2 := {x G n | dist 0O (x, fi c ) > 2M} . (2.13) 

We observe that the £°°-distance of the points in Q x to the boundary d£l of il is 
less or equal to 2M, so fl^ffli Q(2M), where Q(£) = {-I, . . . , £} d + LZ d . 
Hence 

< |dfi|-|Q(2M)| = (4M+ |<9fi| , (2.14) 

and therefore 



= ^p(x) + ^p(x) < (4M + l) d 5|<^| + (2.15) 

rrGSl rrSf2i £6^2 :rgSl2 
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since p < 5 on f2. Eq. (12.151) and Lemma |3~T1 yield 

Tr{K^} > Tr{[P^Pn]-} (2.16) 

- (Ad 5 1 ' 2 + n {am + i) d 5) \an\ - ii J2 p( c 



Next, we apply Lemma l3~2l which asserts 

Tr{[Kp + Up\-} > Tr{[Pn(K, + Up)P n ]-} - ^1^1- (2-17) 

Denoting by x '■= ^[Pn(K^ + Up)PQ < 0] the orthogonal projection onto the sub- 
space of negative eigenvalues of Pq(K^ + Up)Pn and p x (x) := Xx,x its diagonal 
matrix element, we observe that 

Ti{[P Q {K^ + Up)Pn}-} = Ti{Pn{K^ + Up)Pux) (2.18) 

= Ti{P n K,P nX } + UJ2px(*)p(x)- 

By Lemma [331 the density p x is bounded below on fi 2 by the universal constant 
k(ijl) > defined in (13.191) . Therefore 

Tr{[K M +Pp]_} > Tr{[P^P n ]_} - ^ + ^f/p». (2.19) 

Adding up (12.161) and ( 12.191) . we obtain 

^ Z) (7) > Tr{[P^P n ]_} + Tr{[P^X]-} 
- [Ad 6 1/2 + fi (AM + l) d 5 + — } \dn\ 
+ U P( X ) ~ PPW} > (2.20) 

and Lemma l3~6l further yields 



4^(7) - Hl K ^-} > {«(/i)-4^ 1 / 2 - / i(4M + l)^-^}|^| 

+ (2.21) 
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We choose 

. r a(/i) 2 «(//) k(p)\ 



I (12oT) 2 ' 3/i(4M+ l) d ' 2 J ' (2 ' 22) 



and we observe that if 



r 2/i 24 d 2 ) 

U _ U*(p,a(p)) \ 8*(p, a(p)) ' a(//)) J ^ ' ^ 

then our choice for 5 fulfills the requirement 5 > 2pU~ 1 . Moreover, Eqs. (12.221) 
and (12.231) imply that 

4d^ + Ai(4M + l)'*+8* < ^) + ^) + ^) < (2 . 24) 

(7 o 6 6 

We further set Q' 2 := {x G ft 2 | p(x) < ^j} and f^' := {x e Q 2 \ ^ < p(x) < 5}, 
so £l 2 is the disjoint union of Q! 2 and f^', an d by the definition (12.101 ) of p, we have 
that 

{n(p) U p(x) - p p(x)} (2.25) 
> Y^{k{p)U - p) p{x) + £ | {«(//)- 25} > 0, 

since 5 < \n{p) and £/" > 2p/8*{p,a{p)) > p/k{p). Eqs. (12.241) and (12.251) in- 
sure that the right side of (12.211) is nonnegative, which immediately implies Theo- 
rem [Lll 

Theorem ll.2l is obtained by substituting the explicit value of a(p) from (13.601) 
into d2~23l and using L* (//) from (l3~60l) . QED 



3 Auxiliary Lemmas 

In this section we state and prove the lemmas used in the proof of Theorems ll.il 
and ll.2l in Section |2| 



3.1 The Region Q c of High Density 

In this subsection, we estimate Tr{i^ t 7} from below. We are guided by the 
intuition that 7 is essentially localized on fl c . 
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Lemma 3.1. 

Tr{K^} > Tr{[P^Pn]-} - 4d8 l / 2 \dQ\ - /i^p(x). (3.1) 
Proof. Inserting 11 = P n + Pq into Trji^ 7}, we obtain 

Tr{^ 7 } = T^i^P^P^} - 2ReTr{P f }AP n7 } + Tr{K,P a jPa} 
> Tr{[P^Pn]-} - 2 Yl A ^yhy^\ ~ ^{P nl P n } 

= Tr{[P£K„P£}„} - 2 A ^l7^l " /^E^)' 

(3.2) 

where we use that -A > 0, that P^AP^ = P^AP^, and that < 7 < 1. 
The latter also implies that p(y) = 7^ < 1, for all y E A. Thus, if x E d£l 
and y E fi c , the Cauchy-Schwarz inequality yields 17^. | < ■ ^x,x < 

Moreover, if x E d£l, y E fi c , and A x>y ^ 0, then y is a neighbor of x, and we 
obtain 

E A ^v\^\ < 51/2 E E = 2d5 1/2 |0n|, ( 3 - 3 ) 

x£dQ,yen c x£dQ y£A:\x-y\=l 

which completes the proof of (13.11) . QED 

3.2 Decoupling the High and Low Density Regions 

This subsection is devoted to showing that Tr{[i^ M + Up\_} essentially agrees 
with the corresponding eigenvalue sum Tr {[Pq(K^ + £/p)Pn]_} for the operator 
localized on fi, the reason being that tt c is a classically forbidden region since 

-p + Up> \U5 > on Vt c . 

Lemma 3.2. 

Tr{[/^ + f/p]_} > Tr{[P n (^ + [/p)P n ]_} - (3.4) 

Proof. We wish to apply of the Feshbach projection method. To this end, we first 
observe the following quadratic form bound, 

P^(K^ + Up)P£ > P^(Up-jl)P£ > ~U5P£, (3.5) 
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for any p G [0, p], since p > 5 on il c and 5 > 2pU \ Thus, Pq (f^ + Up)Pq is 
positive and invertible on RanP^ , and moreover, we have that 

P n A P£ [PfriKfi + Up)P£] - l P^ AP n < ^ Pan A P^ A P 0Q . (3.6) 
For y E Q c and / 6 C A , the Cauchy-Schwarz inequality implies that 



(f\ Pan A \ A P dn f) = |(APan/)M| 2 = ^ f(x 

x£dQ,\x—y\i=l 

< ( £ |/(aOr) ( £ l) = 2d £ |/(a0| 



(3.7) 



a'£an,|x— j/| i=l a;GA,|:r— y 1 1 =1 x£dQ,\x— y 1 1 =1 

which, by summing over all y G fi c , yields 

(/|PanAP^-AP 9n /) = £ </| P OT A 11, A P dn f) (3.8) 

j/en c 

< m £ • ( £ 0} 

xean yeA,|x-y|i=i 

< 4rf 2 £ l/WI 2 = 4ci 2 (/|Pan/>. 

xean 

(We thank D. Ueltschi for pointing out (I3.7I) - (I3.8I) to us.) We conclude that 

Pn A P^ [p£(Kji + Up)P^] -1 P^ AP n <^ P dn . (3.9) 

The invertibility of P^(Kjx + Up + e)P n L on Ran Pq implies the applicability 
of the Feshbach map, for any e G [0, //]. I.e., for any e G [0, /i], 



P(e) := F Pn [K^ + e + Up\-ePn (3.10) 
= P n (Jf„ + Up)Pn - Pn A J# [P^ + e + Up)P£] _1 jtf A P^ 



is a well-defined matrix on RanP^, and the isospectrality of the Feshbach map 
guarantees that — e G [— //, 0) is a negative eigenvalue of + Up of multiplicity 
m(e) if and only if — e is an (nonlinear) eigenvalue of P(e), i.e., if the kernel of 
P(e) + e, as a subspace of Ran Pn, has dimension m(e). Note that P is monoton- 
ically increasing, as a quadratic form, in e > 0. In particular, 

P(e) > P(0) > Pn(K^ + Up)P n - jfjPixi, (3.11) 
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additionally taking (13.91) into account. 

We claim that, for all A G (0, oo), the number of eigenvalues of + Up 
below —A is smaller than the number of negative eigenvalues of F(X) + A, 

Tr{l[K„ + Up<-X]} < Tr n {ll[F(A) + A<0]}, (3.12) 

where Tr^ denotes the trace on RanP^. Both sides of Eq. (13.121) are zero and 
thus fulfill the claimed inequality, for A > p. Assume that (13.121) is violated, for 
some A G (0, oo), i.e., that A* := inf{A G (0, oo) | Eq. (13.121) holds true} > 0. 
We show that this assumption leads to a contradiction. Obviously, — A* must be 
an eigenvalue of + Up, and hence also of F(X*), of multiplicity m(A*) > 1, 
because only then the left or the right side of ( 13.121 ) changes (increases, in fact). 
Moreover, Eq. ( 13.121 ) holds true for A = A* itself, i.e., the infimum in the definition 
of A* is a minimum. Hence, for all sufficiently small e > 0, the definition of A* 
and the monotony of F(e) in e yield 

Tr{l[K^ + Up< -A*]} < T^{1[ J F(A*) + A* <0]} (3.13) 
Tr{ll[_K" M + Up < —A* + e\] > Tr Q {ll[F(A* - e) + A, - e < 0] } 

> Tr n {J[F(A,) + A, -e<0]}. (3.14) 

Choosing e > so small that — A* is the only eigenvalue of K^ + Up in the interval 
[ — A*, — A* + e], we hence obtain 

m(A*) = TY{1[0< K^ + Up + X* < e]} 

= Tr{l[K^ + Up < -A* + e]} - Tr{l[K^ + Up<-K}} 

> Tr n {l[F(X m ) + A* < e]} - Tr Q {ll[F(A*) + A* < 0]} 

= Tr n {]l[0 < F(A*) + A* < e]} = m(A*) , (3.15) 

arriving at a contradiction, which proves (13.121) . for all A G (0, oo). From (13.121) 
and (13.1 lL we finally conclude 

POO 

Tr{[i^ + f/p]_} = -/ Tr{l[K„ + Up<-\]}d\ 

Jo 

POO 

> - / Tr Q {ll[F(A)+A<0]}. 

Jo 

/oo 
Tr Q {ll[F(0) + A <0]} (3.16) 

= Tr{[F(0)]_} = Tr n {[F(0)]_} 

> Tr{[P n (K, + Up)Pn]-} - ^Tr{P dn } . 
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which is the assertion of Lemma I3T21 QED 

3.3 The Electron Density in the Bulk 

In this subsection we consider the spectral projection 

X := l[Pn(K ll + Up)P n <0] = ll[Pn (-A - p + Up) P n < 0] (3.17) 

of Pq (—A — p + Up) Pq, onto its negative eigenvalues. Writing A^ := Pq A Pq, 
i.e., {Aa) x ,y = A XjV , for x,y G Q, and = 0, otherwise, and V = J^n^x) ' 
1L. := /xP n — UpPn, we have that 

X = 1[-A n - V < 0] and Vx G : -/4 < V(ar) < // , (3.18) 

due to the definition (12.91) of p. Naive semiclassical intuition tells us that, for 
x G fi, the particle density p x (x) : = Xx,x corresponding to the one-particle density 
matrix x should be bounded below by the particle density of the Fermi gas given 
by the one-particle density matrix 11 [—A < p/2]. The purpose of this subsection 
is to prove such a bound (up to a constant factor) where it can be expected to hold, 
namely, for those points x that are sufficiently far away from the boundary of £1. 

Lemma 3.3. Let < p < Ad, define M : = M* := 12(y ) 2 . Suppose that L obeys 
L > 2M and that x G with dist^x, dVl) > 2M. Then 



pT r i -2d 



p x (x) > n(p) := gJX^^ [l + 2 1n(2)(d- 1 + l)+ln(4^- 1 ) 



(3.19) 



Proof. For any (3 > 0, we note that the map R n -»• R, W i-> (e^^ -1 ^)^ is 
monotonically increasing in VT. Namely, as Tq = PqTPq has nonnegative matrix 
elements, so does e eAn , 

OO 

(e^A = e - 2de (e eT ") = e~ 2de Y £ - (T&) >0, (3.20) 

V / w,z V / W,Z / -i h\ V ii ->w,z — ' v 

fe=0 

for all w, z G Vt. So, if n is an integer and W, W G M. n with W^) < W(z), for 
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all z E f2, then we have that 

n 



^ Yl { II ( e/?An/n ), ! z eP&M/A 
2i,...,« n _ien j=i 

= Q e /3An/n e /3W7n]^ ? (3 21) 

V / z , z n 

for all z , z n E £1. Setting z := z n := x E and taking the limit n — > oo, the 
Lie-Trotter product formula and Eq. (13.211) imply that 

/ e -/J(-A n -W0N < / _^(-A n -W)N (3 22) 

V / a;,x — V / x,x 1 v ' 

indeed. In particular, 

e M*( e ^n) < ( e -/»(-An-V)^ (3 _ 23) 

since V > |/ionfi. On the other hand, — A^ — V > — \x and x _L (— A^ — V)x ± > 
0, as quadratic forms. The spectral theorem thus implies that 

xe -p { -A Q -v) x < xe to x = e?»X, (3-24) 

x ± e -/J(-Ao-V) x ± < x ± < . (3 . 25) 

Putting together ( 13.231) . ( 13.241 ). and ( 13.251 ). using that X and — Aq — V commute, 



we arrive at 

V ) x,x — V / x.x 



< e*X» + 1- (3-26) 
Solving for p x (x) = Xx>x , we therefore have 

Px (x) > e-^[(e^) X}X - e-M*\, (3.27) 

for any x E Q and any j3 > 0. 

Next, recall that Q(M) = {— M, . . . , M} d + LZ d = {y E A : |y|oo < M} is 
the box of sidelength 2M + 1 centered at E A. Since dist oc (x, <9f2) > 2M, by 
assumption, we have that 

Q(M) - z + x C fi, (3.28) 
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for all z E Q(M). By Lcmmal3~4l this inclusion implies that 

(exp^Ao])^ > (exp[/3A Q{Myz+x ]) x/x = (exp[/?A Q(M) ])^ , (3.29) 
and by averaging this inequality over z E Q(M), we obtain 

(exppAo])^ > £ (exp[/3A Q(M )])^. (3.30) 



Now, we apply Lemma 1331 and arrive at 

1 e -d/3/M r 

f2e~ l3/M r /2 \ d 
= ( — - — J exp [-4(3 sin 2 (t)] dt J , (3.31) 

where u;(A;) = u(—k) = J^ =1 2{l — cos^)} = J^ =1 4 sin 2 (fcj,/2). Choosing 
/? > 1, we observe that i/^e"* 2 ^ > \ £ e~* 1 dt = ^erffvr] > \. Using 
this and sin 2 (t) < t 2 , we have the following estimate, 

2e -/3/A/ r-K/2 e -P/M i rs/$* 2 e -p/M 

— I «ph/> -»•(«)]*> -^•-y o 

(3.32) 

Inserting this estimate into (13.311) and then the result in (13.301) and (13.271) . we 
obtain, for any /3 > 1, that 



p x (x) > e"^ 2 



-d/3/M 



/\A fid/2 



1 ; ll6ed ri 



(3.33) 

where r := /3/i/<i. Note that if we require r > 4 then /3 = rd/jU > 1, since 
fi < Ad. We may thus replace (3 E [1, oo) by r e [4, oo). Our goal is to choose r 
such that 

P i- \ d/2 

^ ^ > 2 ^ (3.34) 



16ec? r 



r - ln(r) > F := 1 + 2 ln(2) (j + l) + In (jf) . (3.35) 
Note that, due to /i < Ad, 

2.38 < l+2 1n(2) < F < 3 In (16 d ^) . (3.36) 
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We choose r := Y + 2 ln(Y) and observe that Y > 2.38 insures r > 4.11 > 4, as 
required. Moreover, with this choice, we have 

T-ln(r)-Y = 2m(Y) - In [Y + 2In(Y)] 

> ln(F) - In [1 + 21n(F)F- 1 ] (3.37) 

> ln(F) - 21n(F)F- 1 = 21n(Y)(~ - i) > 0, 

using that ln(l + e) < e, for e > 0, and Y > 2.38 > 2. Thus, (13351) and (l3~34h 
hold true. Additionally, we observe that Y < 3 ln(— ) and 



r < F-maxU + 2 

r>0 



{ 1 + 2 (v")} = ( 1 + 2 / e ) y ^ 2F 



(3.38) 



insures that ^f<^f ln(^) < 12(^) 2 < M» < M. This, in turn, yields 



exp 



'4(^2 

2dr 



MfJL 



> 1 



(3.39) 



and by inserting (13.391) and (13.341 ) into (13.331) . we arrive at 



p x (x) > e 



-TlJ 



1 1 



£2d+l e d t 



1 + 2 ln(2) {dr 1 + 1) + In (4rf^ 1 ) 



-2d 



. (3.40) 
QED 



Lemma 3.4. Let A, B C A, with A C B, and denote A A := P4AP4 and A B '■= 
P B AP B . For all x e A and all (3 > 0, 

(exp[/3A A ])^ < (exp[/3A B ])^. (3.41) 

Proof. We first define the nearest-neighbor hopping matrix T on A by T WiZ := 1 
if |iu — 2f|i = 1 and T W)Z := 0, otherwise. For a given subset C C A, the 
matrix Tc ■= PqTPc denotes the hopping matrix restricted to C. Note that 
Ac = Tc — 2dPc is the difference of the two commuting matrices Tc and 2dPc- 
Hence, for x G C, 

(exp[f3A c ]) xx = (exp[f3T c ] exp[-2d/9P ]), x = e~ 2 ^ (exp[(3T c }) x x . 

(3'.42) 

Due to this identity and the fact that x e A C B, Eq. (13.411) is equivalent to 



(ex P [/3T A ])^ < (exp[/3T B ]) . 



(3.43) 
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Now, < { Ta)w,z < (T B )w,z, and hence (TJ)*^ < C^i*)*,*. for a11 intergers n. 
Thus, (13.431) follows from an expansion of the exponentials in Taylor series, 

00 an 00 /3n 

n=0 ' n=0 

(3.44) 
QED 

Lemma 3.5. Le? Q = {— m, . . . , m} a! C Z d be a cube. Denote by Aq the nearest- 
neighbor Laplacian on Q, i.e., A Q = PqAPq = -2dP Q + T Q , T Q := PqTPq, 
andT xy = — y\i = 1). Then, for all (3 > 0, 

— ^(expI/JA,])^ > 0.45) 



where uo{k) := Ylt=i^{^ ~ cos(fcy)}. 

Proof. We may pick an even integer r, choose L := r ■ (2m + 1), and identify Q 
with Q + LZ d C A. (Note that the statement of the lemma makes no reference 
to the Hubbard model analyzed before, and for the purpose of the proof, L can be 
taken an arbitrarily large integer multiple of 2m + 1.) Given s G Z|f, we define 
Q(s) := Q + (2m + l)s and observe that the family {Q(s)} sgZ d of cubes define 
a disjoint partition of A, i.e., 

A = |J Q(s) and Vs ^ s' : Q(s) n Q(s') = . (3.46) 



Hence 

A := J2 A QM (3-47) 

is the sum of translated, but mutually disconnected copies of Aq. We observe that 
Tr{exp[/3A]} (3.48) 
= Z)(exp^]) x>iB = ED ex P^]) 2+ ( 2m+ i)v + (2 m+ i) s 

= SS( eX P^ A «wO, + (2m + l) Sl ,+(2m+l) S = ^ £ ( ^P^qI) *,* ■ 
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As an intermediate result, we thus have 

E ( e M(3^Q}) z>z = Tr{ exp[/3A]} , (3.49) 

since |A| = L d = r d |Q|. 

Next, we translate A by the elements of Q, i.e., for rj E Q,we introduce A^ 
on C A by 

A (,) := A Q(g)+T; = J2 AQ + , +(2m+ i) g • (3.50) 



Of course, A^ is unitarily equivalent to A. We observe that 

l^f E i( "> = E = - M ■ v + E ■ (3.5D 

where, for w, z E A, 

V — f / u;,z — • 

yeA yeA 

= \(Q + w)n(Q + z)\ -T W:Z = 2m(2m+l) d - 1 T w ^ z , 
since T W:Z ^ only if w — z are neighboring lattice sites. Hence, 



A 



_ > A (T/) = -2d • 1 c a + T = 1 c a + 

Q ^ C 2m +1 2m +1 c 2m + 1 

1 1 veQ 

> -— ■ 1 c a + A (3.53) 

m 

where A < is the nearest-neighbor Laplacian on A (with periodic b.c). This 
and the convexity of A i— > Trje' 3 " 4 } therefore imply that 

Tr{ex P [M]} = ^|^Tr{exp[M (T/) ]} > Tr|ex P [A^ A (,)]| 

> e^ d/m Tr{exp[/?A]} . (3.54) 

We diagonalize A by discrete Fourier transformation on C A . The eigenvalues of 
—A are given by u(k), where k E A* = is the variable dual to x E A. Since 
| A* | = L d = \Q \ r d , we therefore have 

1 1 ^ p —/3d/m 

— 5>p[/?A Q U = -Tr{exp[/3A]} > — — £ e HM*) . ( 3 .55) 
zeQ 11 11 fceA* 
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Inequality (13.551) holds for every L = r{2m + 1), and hence also in the limit 
L — > oo. Since the right side of (13.501) is a Riemann sum approximation to the 
integral in (13.45b . this limit yields the asserted estimate (I3.45I) . QED 

3.4 The Discrete Laplacians on Q, £l c , and their Eigenvalue 
Sums 

In this final subsection, we compare the sum of the eigenvalues of 

-A := P n (-A)P n + P^(-A)P£ (3.56) 

below fx to the sum of the eigenvalues of —A below ll, where Vt C A is an arbi- 
trary, but henceforth fixed, subset of A, and Vt c : = A \ is its complement. To 
this end, we introduce the difference of these eigenvalue sums, 

SE(ji,0) := Tr{[-A- M ]_} - Tr{[- A -//]_} (3.57) 
= Tr{(-A-/i)P_} - Tr{(-A-//)P_}, 

where P_ := 11[-A < /j] and P_ := 1[-A < fj]. We further set P + := Pf 
and P + := Pf . Since P_ commutes with Pq, we have that Tr{(— A — ll) P_} = 
Tr{(-A -/i)P_}, and thus 

6E(fi,Q) = Tr{(-A-/x)(P_-P_)} (3.58) 
= Tr{[-A-//]_(P_-1)} + Tr{[-A-/.] + P_} 
= Tr{[A + /U ] + P + } + Tr{[-A- M ] + P_} > 

is manifestly nonnegative. The derivation of a nontrivial lower bound on 5E(ll, f2) 
of the form 5E(ll, f2) > a(Li) \dQ,\, where a(ii) > is a positive constant which 
depends only on ll and the spatial dimension d > 1 (but not on Q), is a task that 
was first addressed by Freericks, Lieb, and Ueltschi in [24|. Shortly thereafter, 
Goldbaum [[231 improved the numerical value for ct(/i) > 0, especially if ll is 
close to 2d. As a consequence of the estimates in [24]|25j, we have the following 
lemma. 

Lemma 3.6 (Freericks, Lieb, and Ueltschi (2002), Goldbaum (2003)). 

(i) Let | < ll < Ad. There is L*(li) < oo and > such that, for all 

L > and all subsets HQ, 



5E(li,Q) > a(n) \dn\ 



(3.59) 
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(ii) Let < fj, < 2, and define 

l^d-ll (2+d)/2 4?rrf 

:= 2W(2vr)^(4rf)^ ^ L * (/i) := ■ (3 - 60) 

where \S | ?5 surface volume of the d-dimensional sphere. Then, for all 
L > (fj,) and all subsets S]CA = Z^, we /zave 

8E(fi,Q) > a(n) \dQ\ . (3.61) 

Proof. We only give the proof of (ii), which amounts to reproducing the proof of 
Lemma 3.1 in ll24l . By {V>fc}fceA* ^ C A we denote the orthonormal basis (ONB) 
of eigenvectors of A, i.e., 

M*) ■= |Ar 1/2 e"^', k e A* = — Z£, (3.62) 

and we have that —Aip k = uu(k)ip k , with u(k) = Ylt=i ^{1 — cos(A; I/ )}. Evaluat- 
ing the traces in Eq. (13.581) by means of this ONB, we obtain 

5E(n,tt) = {[»-u(k)]+(ipk\P+ipk) + [u{k) - fi]+ (ip k \ P-ipk)}- 

> J2^- u ( k )]+^\P+^k)- (3.63) 

fceA* 

Let C C A be an ONB of eigenvectors of A, i.e., — A(pj = ejtpj. For any 

k £ A* and 1 < j < |A|, we observe that 

(e^ - u{k)f\{M^)? = |(^|(A-A)^)| 2 

= \^ k \(P n AP^ + P^AP n )^)\ 2 

= \(P n AP^ k \ Vj )\ 2 + \{P^AP^i)\ 2 

> \{P dn AP^ k \ Vj )\\ (3.64) 

using that either PnV^j = or P^ (fj = and that P n AP^ = P 9 qAP^. Since 
\ej — uj(k)\ < Ad, Eq. ( 13.641 ) implies that 

(4d) 2 |(^->| 2 > \(b k \^)\ 2 , (3.65) 
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where b k : = PaoAP^-^fc is the boundary vector that plays a crucial role in ll24ll . 
By summation over all j corresponding to eigenvalues e 3 - > fi, we obtain 



(ipk\P+ipk) > {Ad)~ 2 (b k \P + b k ) 



(3.66) 



for all fee A*. Next, the convexity of A h-> [A]+ and the fact that P + = 1[— A > 

l4 > (4rf)- 1 [-A-/i] + yield 



(h\P + h) > L(b k \[-A-fi] + b k ) > l^K-A-^fefc)] 



Ad 
1 



(3.67) 



Now, for any x G <9f2 there is, by definition, at least one point x + e E f2 c , with 
|e|i = 1. Since b k is supported in <9f2, we have b k (x + e) = 0, and thus 

(6 fc |(-A-^)6 fc ) = J2 { J2 Mx)-b k (x + e)\ 2 - ii\b k {x)\ 2 ] 

xedtt |e| i=l 

> (1-fj) Mx)\ 2 = (l-n)\\h\\ 2 . (3.68) 

xedn 



Inserting (I3.66I )- (I3.68I ) into (I3.63I ). we arrive at 



(Adf 



^2\pL-u(k)]+ \\b k \ 



(3.69) 



keA* 



Next, we use that in the sum in ( 13.691 ) only those keA* contribute, for which 

u(k) = 2~2t=i 2{l ~~ cos(k u )} < |, as < // < 1. This implies that cos(k u ) > |, 
for all v E {1,2, ... ,d}. Hence, for these k, we have that 



\h\\ 2 



cedfl cr=± v=\ 

d 



> T^J ^ ( J] S COS (M ^ + ^ G ^ 



xedn <j=± v=\ 



1 1 xean 



\dtt\ 
4|A| ' 



(3.70) 
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since there is at least one choice for (a, v) such that x + oe v £ Q c . Inserting this 
estimate into (13. 691) . we obtain 



(3.71) 



(3.72) 



Now define q : T d — > A* by the preimages 

q-\k) := k + 

for A; £ A*. In other words, given £ £ T d , the point £ A* is the closest point 
to £. In particular, |f - g(f)U < f . which implies that - < 

by Taylor's theorem. Hence, 

ixi5> -«<*>]♦ - It"-^®)]^ 



£ if -"(*)]- 



7T 7T 

L'L)' 



> £^-2^'-^)]^. (3.73) 



Since, by assumption, ^ < ^ = | and cj(£) < £ 2 , we have 



L2 ? 



2 d / 2 d(d + 2 



- /' 



Inserting (I3.73I) - (I3.74I ) into (13.711) . we arrive at the asserted estimate. 



l+(d/2) 

(3.74) 
QED 
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